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Abstract 

We derive a non-perturbative equation for the large scale structure power spectrum of long- 
wavelength modes. Thereby, we use an operator product expansion together with relations be¬ 
tween the three-point function and power spectrum in the soft limit. The resulting equation 
encodes the coupling to ultraviolet (UV) modes in two time-dependent coefficients, which may 
be obtained from response functions to (anisotropic) parameters, such as spatial curvature, in 
a modified cosmology. We argue that both depend weakly on fluctuations deep in the UV. As 
a byproduct, this implies that the renormalized leading order coefficient(s) in the effective field 
theory (EET) of large scale structures receive most of their contribution from modes close to the 
non-linear scale. Consequently, the UV dependence found in explicit computations within stan¬ 
dard perturbation theory stems mostly from counter-term(s). We confront a simplified version of 
our non-perturbative equation against existent numerical simulations, and find good agreement 
within the expected uncertainties. Our approach can in principle be used to precisely infer the 
relevance of the leading order EET coefficient (s) using small volume simulations in an ‘anisotropic 
separate universe’ framework. Our results suggest that the importance of these coefficient(s) is 
a ~ 10% effect, and plausibly smaller. 
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1 Introduction and Summary 


The measurement of anisotropies in the cosmic microwave background (CMB) by the Planck 
collaboration has provided invaluable information about the origin of the seed of structure [T], 
plausibly created during an early phase of accelerated expansion [2]. While this is a remarkable 
achievement, many questions still remain open regarding the nature of an inflationary cosmology. 
Perhaps the most outstanding one is whether a slowly-rolling, weakly coupled, fundamental 
scalar held played the role of the Higgs mechanism for the early universe; or if the universe 
chose a different path, such as dynamics at strong coupling or supersymmetry, e.g. [HHS]. These 
possibilities remain viable candidates to play a role in the early universe, and have not been 
considerably hindered by the Planck data [3]. In particular, the existent bounds on (equilateral) 
non-Gaussianity [6] are still above some well-motivated thresholds [3lE]. 

Large-scale structure (LSS) surveys are due to become the next leading probe for precision 
cosmology, providing new information about the history of the universe through ambitious ob¬ 
servational programs, currently under way to make very precise measurements. Therefore, the 
analytic control for the largest number of modes in the process of structure formation will allow 
us to place better constraints on cosmological parameters, and to further our understanding on 
the seed of structure. Unfortunately, quantitative predictions even for the simplest of the cos¬ 
mologies is a daunting task, due to the non-linear nature of dark matter clustering. Moreover, 
even if one adopts a numerical approach, including baryonic matter is far from straightforward. 

Hence, motivated by the arrival of large amounts of data, a renewed interest on analytic 
methods in LSS has resurfaced. These can be broadly classihed in two categories. On the one 
hand, studying the impact of long-wavelength (IR) fluctuations on modes around the scale of 
baryon acoustic oscillations, see e.g. [H1 - I15| and references therein. On the other hand, consider¬ 
able efforts have also been directed towards understanding the impact of short-wavelength (UV) 
perturbations. Their treatment is further complicated by the non-linear dynamics on short scales. 
The effective field theory (EFT) of LSS [l2l[T6l|T7] has emerged as a useful tool to parameterize 
the imprint of UV modes in long-distance observables, e.g. |18H25] . (See also [26] for a somewhat 
different implementation.) Motivated by various results from V-body simulations [23 EH], there 
have also been attempts to reorganize the perturbative expansion, e.g. IHH] . 

In this paper we pursue another road, and study the soft (i.e. long-wavelength) limit of 
the power spectrum non-perturbatively. This allows us to scrutinize the actual impact of UV 
fluctuations on the power spectrum at large scales. As a byproduct, we assess the UV dependence 
of the leading order -renormalized- coefficient (s) in the EFT, which capture the ‘finite size’ effects 
on long-distance observables daiisKiziEi]. 

The dynamics of large scale structures may be cast in a compact form using the den- 
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sity/velocity doublet ijja ~ (^,0), with 9 = —V • v, which obeys the equation [32] 


dyjipaiq, i) + ^ac{'n)il^c{q, v) = lacdiql"^ - k, qji + k) 'il^dql^ - k, r])dd{.q/‘i^ + fc, ■??) • (i-i) 


Here (fi, 7 ) are known matrices that depend on the dynamics of dark matter clustering. (Notice 
we use slightly different conventions for the vertices, see sec. |3.2| for more details. Everywhere 
we use that repeated indices/wavevectors must be summed/integrated over.) Within the EFT 
framework, extra terms must be added in order to account for smoothed regions. The relevance 
of these extra terms on long-distance observables can be read off from the mismatch between 
solutions to ( 1 . 1 ) and data, and/or simulationsj^ 

When the equations in 0 are solved using standard perturbative techniques, the coeffi¬ 
cients in the EFT approach also include counter-terms. This is necessary to cancel would-be 
divergences introducing e.g. a cutoff. (The so-called ‘loop expansion’ produces ill-defined inte¬ 
grals for different initial conditions naiiT!.) Even after the cutoff dependence is removed by 
the counter-terms, the EFT coefficients do not necessarily vanish. The left over (renormalized) 
contribution is physical, and accounts for the finite size corrections to ( 1 . 1 ). 


To obtain the dynamics of the power spectrum, PabiQ^v), we multiply both sides of the 
equation by ipdq^q) and compute the statistical average ^3]. While appealing, the resulting 
equations are not very useful to solve for Pab{q,r]) in closed analytic form, especially because the 
equation for the bispectrum. Babe, on the right-hand side depends on the four-point function, 
and so on and so forth. In fact, an infinite hierarchy of equations between the n-point functions 
unfolds. It is then desirable to have a non-perturbative expression which would only depend on the 
power spectrum itself, and derivatives thereof. Even though we do not expect such relationship 
to hold in general, there are instances in which non-perturbative expressions may be found in 
specific limits, in particular when one of the wavenumbers is taken to be soft, i.e. g —0. In such 
case, we can use a relation for the density and velocity fields that resembles the operator product 
expansion (OPE) in quantum field theory. Keeping only lower order terms. 


i)a{q/‘^ - k,r])'iljdq/‘^ + k,v) (^fab{k, + gabik, n,r})^'^'ilJL{q,v) H -, ( 1 - 2 ) 

where dL{q,9) is the linearized long-wavelength perturbation in the q ^ 0 limit, and 

fab{k,g,v) = f^\k,r]) + f^\k,r])Q2{g) , (1.3) 

gab{k,g,v) = f^\k,v)Qdl^) + ( 1 - 4 ) 


with g, = k ■ q/{kq), and Qdld Legendre polynomials. This relation reflects our prejudice 
that short scale physics is influenced by short modes in a very specific way. One incarnation of 


^In principle, the mismatch is also due to the vorticity, i.e. w — \/ x v, which we ignored in |T1}. The latter 
may be neglected in the g —>■ 0 limit [32], but matters on short scales, fc 7> g. However, these effects will not modify 
our main results. See sec. [Sj 
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this principle are the relations between correlators in the soft limit, e.g. [Hj. Strictly speaking, 


the expression in (1.2) is only accurate to the order we work here. A more detailed derivation is 
given in sec. 


Inserting this expression in (1.1), multiplying by 'ipciQ,'!]) performing the statistical aver¬ 
age we obtain 


driPabiQ^V) — ^acPcbiQiV) ^bcPaciQ^V) 


-\PL{q,ri) 



C22{v) + 


2 1 
1 0 


Ci2{ri) 


(1.5) 


up to terms which are higher order in q. Here PL{q,r]) represents the linear approximation, and 


C22{v) 

cMv) 


= 

= c{l\v) + cil\v) + cil\v) 


with 


a 


c 

a 


S'(.) 

SA) 


= dvr / dA:/2^2 ^(/c,?7) , 


dTT 


-y / dfe fil\k,r]) , 




Stt 


y y d/s fil\k,r]) , 
^ f dk f[l\k,rj) . 


C[firi) = 


( 1 . 6 ) 

(1.7) 

( 1 . 8 ) 

(1.9) 

( 1 . 10 ) 

( 1 . 11 ) 


Notice that the expansion in small q is only viable when the integrals in (1.5) are dominated by 


hard momenta k ^ q. We will see this explicitly in sec. 3.2 


We arrive then to the desired equation, written solely in terms of Pab{k,'il) and a set of 
coefficients which depend on fluctuations on short scales. However, in order to use this expression 
and match these coefficients with data or simulations, it is useful to relate C22{v) C'i 2 (r?) 

to the power spectrum of hard modes. This can be achieved using soft limit relations between 
(re -|- 1)- and re-point correlation functions at equal times. For the case of spherically symmetric 
soft modes, these were discussed in [T3], see also [33H37] . For instance, for the density contrast 
one finds that the bispectrum obeys 


Bniiq,-k - q,k,r]Y 


q^O 


> PLiq,v) 


^ ^,9 Id 

3 dk 3 dr] 


pK=0 

^11 


ik,v) 


( 1 . 12 ) 


in an Einstein de Sitter (EdS) cosmology, and after angular averaging. Here k = K/{a‘^H^), 
with K the curvature of a hypothetical background cosmology, while P^ (fc, rj) is the full non- 
perturbative density power spectrum in the presence of spatial curvature. Similar expressions 
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hold including the velocity held [Hj (for which PL{q,r]) remains unchanged by construction.) 
The generalization for background cosmologies other than EdS is straightforward. Hence, ex¬ 


pressions like (1.12) allow us to connect {rj) to the power spectrum on short scales, albeit 
a locally curved universe. We obtain, for example. 


m 


cilHv) = + 5 ^: {(T22y iv) 


d 


d 




dr] 


Ok 


K=o 


wher^ 




(1.13) 


(1.14) 


The angular averaged expression does not constrain the I = {1,2} multipoles. In order to 
extract this information we need to include directional soft modes. As we shall discuss in the 
next section, these can be also absorbed into a cosmological background, this time with the 
addition of anisotropies. (See for instance [33] for the case of a plane wave perturbation.) The 
hnal expression for the bispectrum in the soft limit is not incredibly illuminating. Because of 
the anisotropies, one observes two different expansion rates, parallel and perpendicular to the 
directional soft mode q. Moreover, there will be a curvature term, Kh , which only enters in the 


parallel direction. An angular average reproduces the expression in (1.12). 

At the end of the day we will obtain an evolution equation for the power spectrum in the soft 
limit in terms of UV parameters which depend on variations of the spectrum of hard modes in a 
hypothetical anisotropic universe. The (rj) are then expressed in terms of the variance of the 
displacements in the new cosmology, which includes a series of geometrical parameters. 




^ “ 3 / ^ 


ai( 


\,k±,r]±,r]\\,Kii 


(1.15) 


(where we added the superscript ‘A’ to emphasize these are dehned for an anisotropic cosmology, 
and A;||,A:_l denote the projections relative to q.) 

It is then straightforward to show, from the above analysis, that only contribute to 

the (^ 22 ( 12 ) ™ the equation ( |1.5[ ) for the power spectrum at large scales. We will argue 

here that these parameters depend weakly on hard modes beyond the non-linear scale. In turn, 
this implies that the hard modes which contribute the most to the power spectrum in the soft 
limit live near the non-linear scale. 

Strong support for this claim follows from Fig. where we plot ^■^^(/cmax) for our universe 


-without curvature- as a function of a short-distance scale fcmax, defined as 


cr 


ab 


(^max) — Q 


Att 
^ JO 


dk Pabik) , 


(1.16) 


^Notice that the expression for aab is dimensionful, since it represents the variance in the ‘displacements’. We 
could instead, following the Lagrangian-space EFT approach (LEFT) [12], introduce the parameters (-ai^{q,'r]) = 
Jg yp-Pabik, rj) and {q, q) = q^ ^rPabik, q). However, we find that keeping the q^ makes the soft limit 
more explicit, and the distinction is unnecessary (between and when g —>■ 0. 
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Figure 1: Moments in (1.16) as a function of k„ 
presented in p 


c, at zero redshift. Results based on 7V-body simulations 


The dependence on short-distance scales is less relevant at earlier times. Note that 
<^S5 = 0'ii,ase = 0 - 12 , CTgg = 022- 


at zero redshift. We notice that cr^ 2 (^max) and cr| 2 (fcmax) does not vary significantly for fcmax 
greater than the non-linear scale Anl — 0.5/i/Mpc, in contrast to iTf^(/i:niax)- 

While we do not expect the variance to be independent of the cosmology for fcmax Anl, 
the curvatnre -and other geometrical anisotropic structure- induced by the soft mode should not 
significantly affect cr^ on short scales. (These new parameters depend on the properties of a 
long-wavelength mode, ^ 1, with |q| <C Anl-) In other words, we expect 

<^ab{k max ) ^ (^ahik max ) (1-17) 


for fc m ax S> Anl, with (T^^(A:inax) dehned similarly to (1.16). Then, the conclusion we draw from 
Fig.j^would also apply to the variance in a modified cosmology. Namely, the o'^ 2 (i 2 )(^™a,x) display 
weak dependence on short-distance modes beyond the non-linear scale. 


In order to translate this conclusion also to the UV parameters in (1.5) one more step is 
required. That is because not only contribute, but also derivatives with respect to e.g. 

the curvature as in (1.13). However, for modes deep into the non-linear regime, we also expect 
the variance to be less sensitive to curvature and other geometrical parameters. For instance -for 
the case of a spherically symmetric soft mode- one can show (n = K/{a^H‘^)) 


d 


^ < Pab{k,r]) 


(1.18) 


for k S> Anl- This can be inferred from numerical simulations which test the so called separate 
universe approach |391I42| (another name for absorbing a spherically symmetric long-wavelength 
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perturbation into a locally curved background.) Since -on physical grounds- we expect this 
behavior to be a general feature, also in the anisotropic case, we then conclude that both, (722 (t?) 
and Ci 2 {r]), are weakly dependent on modes beyond the non-linear scale. As a byproduct, this 
implies a mild dependence on short modes beyond the non-linear scale for the renormalized 
coefficient (s) of the EFT approach at leading order in derivatives [l2l |T6l [T7j , as we shall discuss. 


The conclusion that the renormalized EFT coefficient(s) at leading order in the soft limit do 
not depend strongly on the UV modes does not necessarily mean they have to be small. In order 


to test the validity of (1.1) as a correct description of the dynamics, hence assessing the size 


of the renormalized EFT corrections, we will perform a truncation of (1.5) which is amenable 
for comparison with existent numerical simulations. Our truncation consists on keeping only the 


(722 (t?) contribution to (1.5). As we discuss, the rationale for neglecting Cuiv) is motivated by the 
fact that it arises from higher multipoles (in a suitable basis), and its contribution is suppressed 
by a relative factor ~ 20 in perturbation theory. Moreover, we will be using the approximation 
(see [13] and [36l |37| for details) 




_ 1 

dr, 

0 1 

K=0 ^ \ 

. 0 

+ 1 / 




(1.19) 


ab 


We do not expect this approximation to hold deep into the non-linear scale. As we argued, 
the curvature dependence should be less important the shorter the scales. However, given the 
information in Fig. the approximation in (1.19) may be in principle more reliable, since the 
relevant modes involved in (1.5) do not live too far from the non-linear scale. This enables us 


to use (|1.19|), and replace the derivative with respect to spatial curvature in (1.13) by a time 

( 1 . 20 ) 


derivative, obtaining 


1 


7 


d 


C22{v) = 12 a22{v) + ^^Tr(^22iv) 


dr] 


in an EdS background cosmology, now with AT = 0. 


The resulting (simplified) equation allows us then to make a direct comparison with traditional 
numerical methods with relatively accurate success. We show that the solution to this truncated 
equation reproduces numerical results within ~ 10% accuracy at zero redshift, which is also 
the expected error due to the approximations. The matching becomes more accurate -percent- 
level- at higher redshifts. Our analysis does not take into account the plausible scale-dependence 
(through a renormalization group flow) of the renormalized EFT parameters. However, our results 
suggest that the importance of the leading renormalized coefficient(s) is at most a ~ 10% effect, 
and plausibly smaller. 


In the next sections we derive all these expressions in more detail. We first derive the coordi¬ 
nate transformation to a locally curved anisotropic background in sec. In sec. [^ we use this 
transformation to relate the soft limit of the bispectrum with the variance on short-scales. Later 
on we use these relations to extract the coefficients in the OPE expansion in terms of integrals 


7 










of the variance in a modified cosmology. We then perform a truncation of (1.5), keeping only 


the C22{v) contribution and using the approximation in (1.19). We solve this simplified equation 
and show that it fares relatively well against numerical simulations in sec. We also discuss the 
errors induced by applying (1.19) and ignoring 6 * 12 (ry). We conclude in sec. with a discussion 
on the implications of our findings, and in particular for the EFT of LSS. Appendix [X| contains 
some numerical estimates for the impact of 012 ( 1 ])■ 


2 Directional soft modes: A locally curved anisotropic universe 

In this section we display the coordinate transformation needed to absorb a directional soft mode 
into an anisotropic background. We will study this map up to order (q ■ x)^ in the soft mode, 
which is necessary to extract the Ci 2 (i]) coefficient. We also discuss the impact in the soft limit 
relations between the bispectrum and power spectrum in a modified cosmology. Here we follow 
similar steps as in [T3], but extend the analysis to include a directional mode. We first study 
the mapping to a locally curved anisotropic universe and then apply this transformation to the 
soft limit of the bispectrum. (See also [13] for a somewhat related discussion within the separate 
universe approach.) 


2.1 Newtonian mapping 

We start from a FRW cosmology in global coordinates that includes a long-wavelength density 
perturbation with wavevector q, 

6l{x, t) ~ + c.c. , ( 2 . 1 ) 


such that „ „ 

In the Newtonian approximation, we have 


_5^(t)e*<7- + c.c. 


ds^ = — [1 -|- 2^l(x, + a^{t) [l — 2<hi(a;, t)^^dx'^ . 

We will search for a coordinate transformation of the form 


( 2 . 2 ) 


(2.3) 


t = tA + f(tA,XA) , 

X = XA + 9itA,XA) , (2.4) 

with coordinate dependent functions / and g, such that one can absorb the long-wavelength 
perturbation into a locally curved Bondi-type metric in a series expansion in q. We added the 
subscript ‘A’ to emphasize this corresponds to an anisotropic cosmology. We start with the first 
non-trivial order. 
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2.1.1 Parallel curvature scale factor 


After the transformation in (2.4) is performed, we expect the metric to take the form 


ds^ = — dt\ + a^t^{tA)dx\ y + a\{tA)dx\ j^ 

- ^ (fit'll£C^ II(t^) + K^{tA)x\^ a^{tA)dx\ , (2.5) 

plus 0(a;^) corrections we will discuss momentarily. Notice that, at leading order in the 
distinction between a± and ay is irrelevant in the term involving curvature. For notational 
purposes, we drop the subscript ‘A’ in what follows. 


The extra time-dependent functions, (a_L, oy, Ar_i_, iFy), are the scale factors and curvature 
parameters respectively. These are defined in the perpendicular and parallel directions with 
respect to the soft wavevector q, 


X ■ q 

^ q , 

II 

(2.6) 

x± = a; — ®y . 

(2.7) 

The form of the time-time component of the metric enforces 



(2.8) 

while the vanishing time-space part yields 


a^9 = Vf. 

(2.9) 

Hence we find 

f da 1 

(2.10) 

and 

[da 1 

»=/ 

(2.11) 


We introduced a time-independent constant of integration for g, but this is not required for /. 
In EdS, where oc , 6 l oc a and H oc we have 


and 




2 V$l , 

9 = -o^?T 79 +c, 


3a2Ff2 

and we consider only this simple case in the following. 


( 2 . 12 ) 

(2.13) 


Applying now the transformation (2.4) in the metric (2.3) one finds, in combination with the 


constraints (2.10) and (2.11), 


~ dx^ ( I — — 2H I 


da 1 


a H 


+a‘^ dxidxj{Vigj -F VjQi 


(2.14) 
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Notice that the term 


da 1 ^ 10 ^ 

+ 2H / —^ —$L, 
a H 3 


is constant in time. Next we expand <h in q • *, 


<h(a;, t) ~ - 




2 


1 


1 


-diit) [l+iq-x --{q-x) --{q ■ x) + 


Likewise the integration constant c can be expanded as 

c = cix + C 2 {x ■ q)x + c^x^q + O(x^), 

which provides an additional contribution to the metric of the form 

dxidxjiyiCj + VjCi) = dx^{ci + C 2 X ■ q) 

+ dxidxjXiqj{c 2 + 2C3) 


+ c.c. 


(2.15) 


(2.16) 


(2.17) 


(2.18) 


We notice that the constant and linear term in ds‘^ can be removed using (2.17). The last 
remaining contribution comes from the time-dependent part of g, which leads to 


(q * dx')^ 

dxidxj{S/igj + Vjgi) 3 2- - ■ 


q^ 


(2.19) 


Garnering the pieces together, we have succeeded in transforming the metric in (2.3) into (restor¬ 
ing the labels) 

(q ■ XAf 


ds^ = —dt\ -|- a^dx\ + 2a^6L{t) ^^ ^xa) -, 


which has the form in (2.5), after the identification 


( 2 . 20 ) 


a± 

= «, 

(2.21) 

K± 

= 0, 

(2.22) 

a\\ 

= a (^1 - (5L(t)) , 

(2.23) 

K\\ 

= 5H‘^a^6L{t). 

(2.24) 


The angular averaging for these expression is slightly subtle. The reason is because the above 
quantities are defined with respect to the soft wavevector, q, which would then be averaged over 


angles. It is clear, nonetheless, that after averaging {q^q^/q^ -3 the expression in (2.20) 

leads to the usual result [HEl! 

(2.25) 

(2.26) 


K = 

3 


and 


UK = a ( 1 - -Sl ) , 


as expected. 
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2.1.2 More anisotropy 


We now include an extra order in q, which is required to express the 6*12 (?/) coefficient in terms of 
derivatives of geometrical properties of the power spectrum in a modified cosmology. Notice that 
at this order anisotropies become important, and the local cosmology is not simply described in 
terms of tidal forces in an isotropic FRW universe. 

Like before, the integration constant c can be used to eliminate certain terms. For instance, 
it can be used to transform terms of the form dx^x^ into {dx ■ x)‘^, etc. After the dust settles, 
we find the following expression for the metric 


(q • xaY" 


ds^ = -dt\ + a^dxl + 2a^6L{t y'^ dxl 

2 

5 o. ,..{q-XAf , 2 , o 2e r^\(^'dxAY 




-dx\ + 


q. 

-(q-xa), 


(2.27) 


q^ q^ 

with the addition of two new terms beyond (2.20[). This can be mapped into a metric of the form 


ds'^ = — dt\ + a? dx\ j_ + (l + &|| QXa^ Oy dx\ y 

- ^ (^11 + ^11 QXa,\\) dx\ , (2.28) 

'-( |2.24[ ), together with 

b\\ = 25L{t) , (2.29) 

/Cy = . (2.30) 

O 

Notice ^yOy —)• f^yo^, at leading order in 5l- These extra terms do not survive after angular 
averaging. As we shall argue, fortunately contributions involving C'12 are subleading, such that 
we may concentrate on the isotropic terms. 


with the identification in (2.21 


2.2 The soft limit of the bispectrum 


In the anisotropic background, the fluctuations V'^(fc,r/) depend upon two scale factors and 
the curvature parameters iLy,/Cy. Also the power spectrum, F’^(- • •), depends not only on the 
absolute value of the wavevector k but on k± and fcy separately. Following the procedure outlined 
in we have for the soft limit of the bispectrum 


Bin{q,-k - q,k,r]) Pliq^q) 




d 

+5^F’n(fc±,fcy,^±,7|h^||) 


dn 


+ 


(2.31) 


with Tyy = log ay and Ky = K \\/The right side is understood to be evaluated at q\\ = 
77_l = q, fcy = k± = k and Ky = 0 after the derivatives have been performed. The ellipses include 
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higher order terms in q, which incorporate derivatives with respect to e.g. /C_l in (2.27). These 


are in principle required to extract the coefficients entering in (1.7) for Ci 2 {t]). To avoid a -rather 
cumbersome- expression we do not write them down explicitly, but they can be straightforwardly 


obtained from (2.28). 

On the other hand, only the £ = 0 contribution is needed to extract C22{'f])- The latter can 
be obtained from (1.12), which follows by averaging (2.31) over angles, i.e. m 


Bniiq,-k - q,k,r]Y 


q^O 


> PLiq,v) 


^ ^,9 Id 

3 dk 3 dr] 

+ 1^ PfAk,r])\K^o 


pK=0 

^11 


ik,v) 


3 dK 


(2.32) 


3 A non-perturbative equation for the power spectrum 


In here we derive all the necessary ingredients that enter in the non-perturbative equation (1.5) 
for the power spectrum at large scales. An intuitive way to derive this equation is to use an OPE 


for two fluctuations at nearby points. This will be done in sec. 3.1 An explicit example of a 
product expansion is given by the non-perturbative relations between correlator functions in the 
soft limit, e.g. [TT|. These relations can be used to fix the coefficients of the OPE in terms of the 


power spectrum on short scales, albeit in a modified cosmology. This will be discussed in sec. 3.2 


3.1 Operator product expansion 

For the theory of structure formation we may think of the gravitational potential, ^{x,r]), as 
the relevant degree of freedom from which a solution for 'ijjaix,'!]) follows, e.g. by deriving the 
displacement fields that follow from the induced forces, i.e. —V<h(a;,r/), as a function of time. 
An OPE is then an expansion for product of V’a fields at nearby points, 

Ya{x,r])Yb{y,ri) + y),r]^ , (3.1) 

O 

as a function of (composite) operators that are built in terms of ^{x,rj) and its derivatives. In 
Fourier space it reads 

i>a{q/2 - k, r])ipb{q/2 + k, y) ^ f^Ak, v) 0{q, rj). (3.2) 

o 

For our purposes it turns out to be sufficient to keep only terms linear in $, because of the 
smallness of the density perturbations (and gradients) in the g —)• 0 limit. However, let us briefly 
mention some caveats of the OPE expansion in our setting. Fist of all, only a handful of examples 
are known to obey an OPE beyond perturbation theory. Moreover, for the theory of structure 
formation, an extra complication arises. This is due to the statistical properties of the initial 
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state. In other words, the existence of different possible realization introduces stochastic terms 
which are not necessarily proportional to products of long-wavelength fieldsj^ However, for our 
purposes, we may ignore these terms, which are known to enter at higher orders in the evolution 
equation for the power spectrum, i.e. jUj. (The relevance of these terms, as one approaches 
the non-linear scale, has been recently emphasized in 12Q1EI].) Furthermore, the operators in 
the right-hand side of (3.2) may themselves involve products of fields at the same point, which 
in turn need to be regularized/renormalized |12j . However, the net effect for these extra terms is 
sub-leading, since they are suppressed by extra powers of Pl- 


Due to the equivalence principle, we know the expansion starts with -at least- two derivatives 
of the gravitational potential. Thus, working in Fourier space and using isotropy, parity, and the 
equivalence principle, we need 


Oo{q,ri) 

= q^^q,q) (x5{q,q) , 

(3.3) 

Oi{q,q) 

= g2g*$(q,rj) oc g*5(q,r?) , 

(3.4) 


II 

(3.5) 

(q, q) 


(3.6) 


where TF stands for trace-free. (The traces renormalize the coefficients of tensor operators with 
fewer indices, and we do not include them for simplicity of notation.) Moreover, for the functions 
of the hard modes we can, in addition to /^°(fe, ??), use a similar splitting into scalar functions: 




oP 

C>tf 




(3.7) 

(3.8) 

(3.9) 


When the dust settles, after performing all the contractions, we wound up with 


V'a(q/2 - k,r])ijjb{q/2 + k,q) (^fab{k, fj,,r]) + Qabik, 


+ 


(3.10) 


where fab{k, /x, r]) and gab{k, /x, rj) are polynomials in [i, up to order with coefficients given by 
a combination of the scalar functions defined above. Higher orders in g arise when expanding 
in higher orders in q. These functions can be re-written in terms of Legendre polynomials as in 


(1.3)-(1.4): 


fab{k,n,r]) = fil\k,r]) + f^\k,r])Q2{g) , (3.11) 

9ab{k,ti,v) = fil\k,r])Qi{n) + fj^l\k,r])Q3{iJ,) . (3.12) 


^These are ‘contact-terms’ in the jargon of quantum field theory. 
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The expression in (|3.10) implies, 


Biab{-q,q/2 - k,q/2 + k) (^fab{k, n,v) + gab{k, PL{q,v) 


+ 


(3.13) 


for the bispectrum in the soft limit, and also with Biab —>• B 2 ab- (Recall is defined such that 
Bu{q,il) = P 22 {q,q) = PL{q,q) at linear order [33].) Notice, as expected, we do not have a 
term proportional to 1/q. This is related to the so-called consistency condition for large scale 
structure correlation functions, which vanishes when evaluated at equal times. This is nothing 
but a reflection of the equivalence principle, see [l5] and references therein. 


3.2 The fluid equations in the soft limit 

The dynamics in 0 leads to 

drj Pabi.q^q) — ^acPcb^Q^q) klbcPac^Q^q) 

+ lbcd{-ql‘i‘ -k,k + ql2) Bacdiq, -ql‘i‘ -k,k- q/2, g) 
+ lacd{ql‘^ -k,k + q/2) BMci-q, q /2 -k,k + q/2, rf) , 

for the evolution equation for the power spectrum [32]. Moreover, 

li2i{k,p) =7ii2(p,fc) = a{k,p) , 
l222{k,p) = p{k,p) , 

with 


(3.14) 


a{k,p) = 

The Qab matrix is given by. 


(k + p) ■ k 

2P 


(3{k,p) = 


O 

2 ^ ‘m 


-1 
1 _L 

dr 


{k + p) k ■ p 

2fcV 


0 -1 


(3.15) 


(3.16) 


(3.17) 


The last equality applies in an EdS background. It is straightforward to generalize idab to other 
cosmologies, such as ACDM [33|, or to include curvature effects [T3|. In our numerical analysis, 
we use the EdS values for the rescaled variables ^|J = (<5, —O/Pf), see e.g. [32]. The error induced 
by this approximation is subdominant. 

The main observation now is the following. In the limit g —?■ 0, the integral in (3.14) is 


dominated by momenta k ^ q. In particular, for a physical power spectrum scaling as P{k) ~ k 
in the soft regime, we have, when k < q, 

d^kP{k)^q^. (3.18) 

This is subdominant, in the g-soft limit, with respect to the contributions from modes with k ^ q, 


scaling as q'^cr^, like we find in (1.5). Moreover, as we see below, is dominated by modes near 
the non-linear scale. 
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In what follows we split the solution to the equation (3.14) in three sets, 


Pab{q, v) = Pab'^iq, v) + Pabiq, v) + Pab^q, v)- (3.19) 

The hrst contribution corresponds to the homogeneous solution, which is identical to the linear 
power spectrum, UaUbPiiq-, v) with Ua = (1,1), when ignoring the decaying mode. The other two 
contributions are sourced by the a- and /3-terms, respectively. We discuss first. 


3.2.1 /3-terms 

The /3 contribution is of the form (the subscript a G {1,2} is arbitrary) 

/ cP’k 

/3{-k + q/2, k + q/2) Ba 22 {-q, ql‘2^ - k, q/2 + k ). (3.20) 

This is quite fortunate, since the vertex function scales like /3 oc q‘^/k'^ for q k, and thus the 
integral amounts to an angular average of the bispectrum, 

-2 r (Pk Ba 22 {-q,q /2 - k,k + q/2) 


2 


( 2 vr) 




(3.21) 


This means, using the OPE, only f 22 \k,q) contributes to the above integral at leading order 
in q. At the end of the day we obtain 


Pabi^^ v) + ^acP^iq, q) + ^bcP^iq, q) = \ Pdq^ q) 


0 1 
1 2 


C22{q), 


with 


C22{q) =4tt J dk f^l\k,q), 

for the contribution to the dynamics from /3-terms. 


(3.22) 


(3.23) 


We could in principle directly use (3.13) to extract f^{k,q) from the bispectrum in the 
soft limit. (After all, this is similar to extracting the (local non-Gaussianity) parameter in 
the squeezed limit for the CMB counter-part.) However, because this case involves the angular- 
averaged bispectrum, we can use the soft limit relations discussed in sec. Using (3.13) and 
(2.31), and performing the angular integration, we recover (1.13) (see also ( |2.32 ) and |14j.l 
Moreover, using the approximation (9 k— |(^ 5 ? + 1)V’^^*^ [E]) we obtain 


(3.24) 


^ ^ 12 2 , X 13 5 2 / ^ 

C22{q) ^ yO-22(^) + y ^^22(^) • 


We will use this result to compare against numerical simulations. 
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3.2.2 a-terms 


While the a-terms appear to contribute at order q, 


q k q‘^ {q ■ 


a(fc + q/2, 9/2 - ^ - 




(3.25) 


since the integrand does not change if both k and q change sign, the leading contribution will 
ultimately scale as q'^, similarly to the /3-terms. Then, using (3.13) we get 


dvPabiq^v) + ^acPS,{q,q) + ^bcPadq,q) = -\ PLiq^ 9) f ^ M ^12(7?), 


with 


and 


Cu{v) = c[l\v) + C[d{q) + C[d{v) , 


(1)/ 


.(2)/ 


(3.26) 


(3.27) 


C[l\q) = 

Att 

"T j 

[ dk f[f{k,q) , 

(3.28) 

cSiv) = 

Stt 

1 dk f[l\k,q) , 

(3.29) 

cilKv) = 

167r j 

J 

dk fil\k,r]) . 

(3.30) 


As we discussed in sec. we can then resort to the locally curved anisotropic universe to 
write Ci 2 {q) in terms of the variance of the displacements, cr(^, and derivatives with respect 
to geometrical quantities. This allowed us to infer that Ci 2 {q) does not depend significantly 
on modes beyond the non-linear scale. However, it requires using the Newtonian mapping for 
a directional mode to O(g^). This is rather cumbersome, as we see explicitly in sec. 


2.1.2 


Moreover, unless simulations can be performed in anisotropic backgrounds, it would be unclear 
how to replace derivatives with respect to these cosmological parameters by information in our 


universe which could be used as input in (1.5). This would be desirable in order to most accurately 


confront our equation with observations, or numerical methods, and asses the importance of the 
EFT coefficient(s) non-perturbatively. There is however another, perhaps more suggestive, way 
to parameterize the a-terms. By a judicious shift in the integration variables we can re-write 
their contribution as follows: 


d^k 

(27r)^ 
d^k q ■ k 


a{k + q/2, q/2 - k) Bai 2 iq, ql‘2 -k,k + q/2) = (3.31) 

[Bai 2 {-q, -ql‘2- -k,k- q/2) - H„i2(q, q/2 -k,k + q/2)] . 


J {2tt)^ 2k^ 

This way we see there is no angular-averaged contribution. Moreover, using an OPE for the term 
between brackets it follows that only /[^^(A:,?/) enters in Cuiq)- (We used a tilde to distinguish 
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from the direct OPE decomposition of the bispectrum.) This is nothing but a rewriting of the 
previous functions. However, because of the explicit angular dependence, it may suggest that 
the a-terms are less important compared to /3-terms in the soft limit. Surprisingly, this turns 
out to be the case in standard perturbation theory. The a-terms are roughly a factor of 20 
smaller than the leading contribution from /3-terms. In the non-perturbative regime, there is no 
justihcation for neglecting these terms, which may still lead to sizable errors. However, since 
most of the modes that contribute to 6*12 (?/) come from near the non-linear scale, we may expect 
only small deviations from the a-terms. Hence, in what follows we work under this hypothesis, 
i.e. P^fj{k,ri) <C P^^{k,ri). As we shall see, this approximation fares well against numerical 
simulations. We will return to this point in sec. [^and present some numerical estimates, for the 
error induced by neglecting the a-terms, in appendix [A} 


4 Numerical results 


We will now concentrate on the non-perturbative equation for the power spectrum in the soft limit 


that follows from the /3-terms, together with (1.19). In other words, we approximate PabiQ^v) — 
Piiq, v) + Pabi^, V), such that 

driPabiQ-iV) ~ ^acPcbiQi V) ^bcPaciQi 

-q^PL{q,r]) f J 2 ) ■ (4-1) 


Furthermore, we parameterize the power spectrum in the soft limit by 

Pab{q, v) PUq, r]) (l + Cab{v)q‘^) , 


(4.2) 


and, after subtracting the linear contribution, the above equation becomes a differential equation 
for the coefficients Cab{r])- 

A simple exercise is to use the linear variance (t |2 ^ right-hand side, 

which yields at leading-order 


^ ^ ^Liv) [ -38 -76 

63 i -76 -114 


- crfinir]) 


-0.60 -1.21 

-1.21 -1.81 


(4.3) 


This must be compared to the direct standard perturbation theory one-loop computation, which 
gives 

( -61 -125 ' 

105 I -125 -189 




- ^Liv) 


-0.58 -1.19 
-1.19 -1.80 


(4.4) 


At this level, the small difference quantifies the influence of the a-terms in the flow equation, 


we alluded before, as well as the error in the approximation in (1.19), which we argued in 
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Figure 2: Numerical results for for the coefficients cn = css,ci 2 = cse,C 22 = cgg as a function of redshift, 
obtained from the (truncated) non-perturbative equation (4.1), using velocity correlation spectra from 
7V-body simulations performed in |38j as input. The black squares in the plot correspond to numerical 
results for the css coefficient for the density power spectrum in the soft limit [22l [23] . The error bands are 
a rough estimate -at the ^ 10% level- of the errors in the non-perturbative equation, induced by ignoring 
Qf-terms and using the approximation in (1.19). The agreement is remarkably good, especially at higher 
redshifts. 


is small within perturbation theory. For the non-linear regime we estimated the error to be 
at the ~ 10% level, which is consistent with Fig. (See also the next section.) The input 
from UV modes in the evolution equation enters through <T22(femax)) the variance of the velocity 
divergence auto correlation up to a given cutoff. To extract this variance we use results from 
Al-body simulations based on a phase-space projection technique [38]. The result is given in 
Fig-i also shown in the lower panels in Fig. For the numerical results we used ACDM model 
with Qm = 0.276, Oa = 0.724, Qh = 0.045, h = 0.703, ug = 0.811, Ug = 0.96. The value of the 
Cabiv) coefficients as a function of redshift is plotted in Fig. For the density coefficients, the 
agreement is remarkably good within the error bars. Unfortunately, we do not currently have 
knowledge of data we could use to compare our predictions for the C 22 and C12 coefficients. 

In order to explore the cutoff dependence in the UV parameters, in the top panels in Fig. 
we show the numerical and perturbative results obtained using two different cutoff scales kmax = 
1(10) h/Mpc. While the perturbative computation depends strongly on the cutoff at the two- and 
three-loop level, the full solution is essentially robust. For the non-perturbative solution, the Cab 
coefficients do not vary significantly when the cutoff changes by an order of magnitude, and the 
main contribution comes from modes with k < 1/i/Mpc. The apparent strong(er) dependence 
on small-scale fluctuations indicated in perturbative calculations, displayed in the standard two- 
and three-loop results, is an artifact. This is also exemplified in the lower panels of Fig. where 
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^max—Oh/MpC 





Figure 3: On the upper panels, solid lines reproduce the results from Fig. Dashed lines show the one- 
loop, dotted lines the two-loop, and dot-dashed lines the three-loop standard perturbation theory results. 
The two figures are given for two different cutoff scales k^ax = 1(10) /i/Mpc, respectively. The solid lines 
are affected at the ~ 2% level by the change in cutoff. For the lower panels, on the left we plot cr^f,(fci„ax) 
(in dashed lines) computed at one-loop as a function of cutoff. On the right, we plot (in dotted lines) 
the two-loop calculation. We notice, once again, the relatively strong dependence on fcmax beyond the 
non-linear scale. This is the case also for cri 2 (fcmax)) especially at two-loops, and also for (T| 2 (fcmax) to a 
lesser degree. On the other hand, the non-perturbative results remain essentially unaltered. 

we see the dependence of the variance cT^^(femax) on fcmax, within standard one- and two-loop 
perturbation theory, contrasted with the non-perturbative results. 

5 Discussion 

We derived a non-perturbative relation for the power spectrum in the limit of long wavelength 
perturbations g —)■ 0, 


dri Pab{q, v) = -^acPcbiq, V) “ ^bcPaciq, ??) 


(5.1) 


19 
































with coefficients C 22 (i 2 ){v) which depend on the variance of velocity-velocity and velocity-density 
correlations on short scales, albeit in a spatially curved anisotropic universe. Let us emphasize, 


once again, that the density-density correlation on short scales does not enter in (5.1), due to 
the structure of the non-linear couplings 7afec- The main ingredients are the OPE together with 
relations in the soft limit between the three-point function and the power spectrum in a modified 
cosmology. 

The presented derivation has important ramifications. Perhaps the most relevant concerns the 
UV dependence of the power spectrum in the soft limit, and the implications for the EFT of LSS. 
At leading order in the soft wavenumber q the EFT of LSS introduces a new parameter, which 
scales as q‘^PL{q, 7) when g —?■ 0. This coefficient includes a series of response functions in LEFT, 
discussed in [12], or may be interpreted as a ‘sound speed’, c^, in the Eulerian approach m- As it 
was originally emphasized in naiiKiT], the effects of the EFT coefficient (s) can be read off from 
the discrepancy between observations, Pohs{k,r]), and the Pfiow{k,T]) which results from solving 


(1.1) without the addition of finite size effects. Alternatively, we can use our non-perturbative 


equation in (5.1) to find the solution for the power spectrum. In order to assess the dependence 
of Z^en oil tho short-distance modes, we introduced a cutoff scale, fcmax- This is the maximum 
wavenumber contributing to the variance of the displacements in density and velocity fields. The 
dependence on /cmax in (5.1) enters through the coefficients 6*22(12) (7)) which are determined 
as a function of the variance (T^^j^2)(^max) •'') in a (hypothetical) locally curved anisotropic 
background cosmology. Hence, for a given /cmax; we can define a (renormalized) parameter which 
tracks the dependence on short scales as follows: 

AP(g, femax 


henikmax) — 


g <C A 


NL, 


(5.2) 


q'^PLiq) 

at a given time. The observed Z^en is given by taking the cutoff to infinity, or in practice 
fcmax S> AnL) i-e. much bigger than the non-linear scale. The dependence on Zcmax then reflects 
the sensitivity to short-distance modes, whereas the importance of these coefficients depends on 
AP{k) = Pohs{k) — Pfiow{k, kraax)- Namely, the mismatch between observations (and/or simula¬ 


tions) and solutions to (1.1). 

In principle the mismatch, AP{k), also measures the effect of vorticity, w = V x v, which was 
ignored in (5.1). The vorticity enters in the density power spectrum through new vertices in (1.1), 


e.g. ‘yscue, and in turn these contribute to the coefficients in (5.1) through terms proportional to 
6'§- (^uoj{kraax)- However, Pu)u}{k) is highly suppressed on long-distance scales, and moreover, 
on short scales it is at most of the same order as the velocity power spectrum [38]. Hence, we 
conclude cr‘^Ukraax) cr‘^ 2 ikmax)■ 

The a'^fj{kmax) are plotted in Fig. as a function of the cutoff, in our universe. Notably, 
cj 22 (fcmax) and (7i2{kraax) show a Saturation on short scales. On the other hand, erf^(Zcmax) ex¬ 
hibits a different behavior. As we argued, the qualitative behavior on very short scales should 
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not depend strongly on details of the background cosmology, which are only relevant at long dis¬ 
tances. For example -if we concentrate on the curvature dependence for a spherically symmetric 
perturbation- we expect 

'^ab(^inax) — Crafe(^max) h)) (^-3) 

beyond the non-linear scale. Moreover, we also expect the power spectrum on short scales to 
obey {k = K/{a?‘H‘^)) 




< Pabik,r]), 


(5.4) 


K=0 


for k ^ Anl- This can be shown to be the case from numerical simulations |391442j . For instance, 
in [l2] the derivative with respect to curvature was calculated, and used to infer the response 
to a long-wavelength perturbation using the separate universe approach. We reproduce their 
results in Fig. for the case of the density power spectrum. As we see in the plot, the numerical 
result lies below the red solid curve, which is constructed using the approximation in (1.19), e.g. 


Q The error is of order ~ 10% around k ~ 0.5/i/Mpc at redshift zero, and the 
separation grows as we increase k. While the results in |42] did not include velocity correlators, 
we expect similar results. On the other hand, we plot 9)7cr^fe(A;max) as a function of femax in Fig.[^ 
We notice, first of all, the weak dependence on hard modes beyond the non-linear scale for the 
22- and 12-component and, moreover, the approximation 9,, log (T22(]^2) (^max) — 2 is relatively 
accurate for fcmax ^ Anl- From these observations we conclude that (5.4) is fulfilled. Therefore, 


the weak sensitivity to hard modes beyond the non-linear scale also applies for the derivative of 
^ 22 {i 2 ) ^ith respect to spatial curvature, which is what enters in (5.1), e.g. (1.13). On physical 


grounds we expect a similar behavior in an anisotropic universe, and the contributions to (1.7), for 


both the variance and the variation with respect to the extra geometrical parameters introduced 
in sec. 12 .1.21 

From the above reasoning, we concluded that the 022(12) (^) vary weakly with scales beyond 
the non-linear regime. This implies 


Tflow(Q') ^max ) — -^Nl) ; ^max > Anl , 


and in turn, for the parameter introduced in (5.2), 


;2 ^ A 'i AP(g,ANL) ^ . 

^renV^max ^ -^Nl) — 9 o 7 ^ ^ -^NL 

rPLiq) 


Strictly speaking the authors in [42] evaluated numerically the response function: 

Gi{k,v) = -\dv'^ogP^p°{k,ri) + | logPn (fc, 77 )! 

o o \k=o 


(5.7) 


(5.8) 


(5.5) 


for the density power spectrum at different redshifts. This is the combination which enters in (2.321. In Fig. we 
compare the results reported in with the (VKPR) approximation implicitly used in |36l I37 |. see also |14| . 

4, 


log Pii (fc, 77 ) 


7 


d,,logPi^ °ik,ri) . 


(5.6) 
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z=0 


z=2 



Figure 4: The response of the power spectrum to a long-wavelength perturbation, Gi(fc, z), at redshifts 
z = 0 (left) and z = 2 (right), as introduced in [42]. The red solid line uses the approximation in (5.6) 
[551157] (see footnote . The results from the A^-body simulations are reproduced here from |42j by 
the black squares. The other curves represent perturbative computations and the halo model. We see 
relatively good agreement between the approximation in (5.6) and iV-body results at low fc’s, as it was 
also pointed out in [T3|. However, (5.6) fails for high fc’s at zero redshift, with an error of order ^ 10% 
near the non-linear scale. The disagreement becomes more dramatic as we increase k, with the numerical 
result much below the red line for k ^ Anl- 



Figure 5: Derivatives of the moments introduced in (1.16) as a function of kmax and at zero redshift. Results 
based on Wbody simulations presented in [38]. Notice that drj log 0 - 22 ( 12 ) (^max) — ^ ^max ^ Anl- (See 
Fig.g) 


Hence, the leading order renormalized coefficient(s) in the EFT frame-work themselves may 
be obtained with information from modes up-to/near the non-linear scale. 

This result is not surprising, and resonates with the statement that virialized scales decouple 
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on large scale dynamics. This was shown in [T^ for the back-reaction on the evolution of the uni¬ 
verse, since short scales contribute through the expectation value of their quadrupole moment |12j . 
which is suppressed relative to non-virialized scales |44) . Likewise, our analysis demonstrates the 
same occurs for the response of this quadrupole to a long-wavelength perturbation, which enters 
at order Our results are then compatible with various numerical investigations of the impact 
of UV modes [271 EH], as well as general expectations from analytical arguments, e.g. miiigiizi. 


This is, on the other hand, in stark contrast with what occurs when standard Eulerian (or 
Lagrangian) perturbation theory is used to solve for (1.1), where Tflow(^,^) is calculated as an 
expansion in the linear power spectrum, involving integrals over all scales. Within a perturbative 
scheme, the EFT parameters are split into a counter-term and a renormalized piece. We find 
that any dependence of the leading order EFT coefficient(s) on short-distance modes beyond the 
non-linear scales develops mainly from the counter-term(s). In other words, the UV sensitivity 
is nothing but an artifact due to the inapplicability of perturbative techniques beyond the non¬ 
linear scale, as opposite to intrinsic UV dependence of the -physical- renormalized parameters 
encoding finite size effects. We see this in detail in sec. where we observe in Fig. [^the cutoff 
dependence in loop computations, in contrast with the non-perturbative results. 


The fact that modes beyond the non-linear scale do not have a major influence on the renor¬ 


malized parameter(s) does not mean that the effect of the new term in the dynamics of (5.1), 
Piiq), is necessarily small. In fact, its relevance relies on how well the solution to (1.1) fares 
against observations. Moreover, on physical grounds, we know that tidal effects are expected to 
contribute on long scales Ha ED. In order to estimate the size of the contribution from the 
renormalized EFT parameter(s) we used a truncated version of our non-perturbative equation. 
Our truncation sets Ci 2 {r]) = 0 on the right-hand-side of (5.1), and uses the approximation 
— ^{dri + 1)V'2'^° applied to (1.13), which implies 


C22{ri) = ^ ( 12 


7 


dr] 


(5.9) 


It is debatable to what extent this expression is valid in the non-perturbative regime. As we 


showed, due to specific features of the equations in (1.1), it seems to be remarkably accurate 
-percent level- in perturbative computations [T3|. However, confronted with V-body simulations 
for the standard ACDM cosmology the accuracy degrades on shorter scales to ~ 10% at z = 0 
and k ~ 0.5/i/Mpc, see Fig. [39H12] . This is indeed expected, as we argued the derivatives 
with respect to curvature and time should not be interchangeable beyond the non-linear scale, 
see (5.4). However, since the modes which contribute to C22{v) th® most live near the non-linear 


^Notice that for virialized objects hvir — 10^. However, using that the potential is roughly constant (of order 
4> ~ 10“®) on all scales, this corresponds to fcvir — IOAnl. Furthermore, the power spectrum turns over near 
k ~ fceq < Anl, namely the horizon scale at matter/radiation equality. This introduces an extra suppressing 
factor, and is the reason we see the turn in Fig. [l] closer to Anl — 0.5 /i/Mpc. 
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scale, the approximation in (5.9) may be more reliable. We see this explicitly in Fig. 


In addition, we also have the a-terms encoded in 012 ( 1 ]), which we neglected in our simplified 
equation. These are subleading within perturbation theory -a factor of 20 smaller- and their 
relevance is unclear in the non-perturbative regime. Since, as we argued, only modes near fcmax — 
Anl contribute to these coefficients, we do not expect the a-terms to modify our results radically. 
Hence, we estimated overall an error in our truncation roughly of order ~ 10% at z = 0. This level 
of precision is consistent with the results presented in sec. see Fig. and also with a simple 
estimate of the impact of Cuir]) in appendix [A} For example, our non-perturbative (restricted) 
solution for the density power spectrum, given by (in redshift space) 


Pssiq, z) = (1 + cs5iz)q'^)PL{q, z) , (5.10) 

with 

cssiz = 0) ~ -21 (Mpc/h)^ , (5.11) 

is within the error of the numerical value [22L 

= 0) ~ -23 (Mpc/h)2 . (5.12) 


The accuracy improves the higher the redshifts. Because of the series of approximations, it is not 
possible to directly extract the size of the renormalized EFT coefficient(s) in the soft limit. We 
can, nonetheless, conclude that their importance in correcting 0 is at most of the same order 
as our error bands, at the ~ 10% level. 


The overall level of accuracy could be significantly improved by fitting the time dependence 
of the response function to a long-wavelength perturbation, Gi(k,z) in Fig. integrated over 
wavenumbers. This is ultimately what contributes to 022 ( 1 ])■ Because of the missing Oi 2 (i]) 
term in our truncated equation, this was not justified at this stage. (In fact, the approximation 
in (1.19) proved to be relatively accurate to compute 022 ( 1 ]), with an error which is similar to 
the one induced from ignoring a-terms.) However, in principle we could also incorporate Oi 2 (i]), 
provided response functions to directional long-wavelength perturbations could be obtained from 
numerical simulations, using an anisotropic separate universe approach. Contrary to a direct 


extraction of the css(z) coefficient in (5.10) from Wbody simulations, where one needs large 
volumes to beat variance in the soft limit, e.g. [221123] . computing the response functions require 
modest simulation volumes [39H42| . This suggest a hybrid analytic/numerical approach to model 
the power spectrum. The numerical input would come from fitting the time behavior of integrated 
response functions, unlike extracting the low-g behavior of the power spectrum. In principle, this 
would allow us to precisely determine the size of the renormalized leading order coefficient(s) in 
the EFT of LSSHThis is left for future work. 


®There is a technical point worth mentioning. The renormalized parameters in the EFT are not necessarily 
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Impact of Ci2 



Figure 6: Estimate for the impact of (^12 ( 77 ) on the coefficients Cab in 
Mill M 2 I £ 3. 


(4.2). 


For definiteness, we used 
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A Estimating the impact of Cu 


In the following we present an estimate for the impact of the coefficient Ci 2 {r]) in (5.1), through 
the contribution from to the power spectrum introduced in sec. 3.2 As was discussed 


constant for all values of q, but may also depend on scale. The origin of this inherited g-dependence is the 
existence of non-analyticities such as logarithms, e.g. g^(logg)". This means, in practice, that the renormalized 
parameters at different scales may differ by large logarithms which need to be resumed. This is achieved through 
the renormalization group flow. On the one hand, our conclusion will not change when we vary q towards the non¬ 
linear scale: For q ~ Anl the renormalized parameter(s) will not depend strongly on short-distance modes beyond 
the non-linear scale. On the other hand, there could still be an important flow between modes with q <C Anl 
and the non-linear scale. We do not expect this scale dependence to have a large impact, but it is nonetheless an 
important factor when hunting for percent-level accuracy. See [SI EH for a somewhat related discussion. 
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in sec. the Ci 2 {ri) coefficient depends on the response of the correlation between velocity and 
density to various geometric parameters, including background curvature, etc. These response 
functions could in principle be extracted from suitable A^-body simulations. In absence of such 


data, and inspired by (5.9), we then explore the following ansatz 

d 

Cuiv) = di o-iaC??) + d2 , 


dr] 


(A.l) 


with some adjustable free parameters {di,d 2 )- In order to estimate their plausible values, we 
evaluate (5.1) within perturbation theory. At leading order, —>■ oc Requiring that 

the left-hand side of (5.1) reproduces the one-loop correction to the power spectrum, and using 


(5.9), one finds 


di + 2d2 = — 7 ^ ■ 
6b 


(A.2) 


Notice this is a significant suppression which, compared to (4.1), precisely accounts for the %-level 


difference between (4.4) and (4.3). 


We can then estimate the impact of Ci 2 {r]) by varying di and ^2 while taking (A.2) into 
account. The results are shown in Fig.[^ The dependence is relatively mild, at the level of 

~ 15% (5%), Ac50~5%(1%), Acee<l%, (A.3) 

at z = 0 (0.5) respectively. 

Notice that, for the cri 2 i'n) extracted from the iV-body simulations, one finds dr^ In ( 7^2 — 2 at all 
redshifts. Hence, the influence of the Ci 2 {r]) term, and related uncertainties, are then effectively 
suppressed due to the small overall coefficient in (A.2). We also find that velocity components are 


even less affected by 6*12 (r/), which can be understood from the structure of (5.1). Note that the 
results presented in Fig.[^ which are obtained neglecting C'i 2 (??), are also within the uncertainty 
bands shown in Fig.j^ In addition, one can adjust the remaining free parameter in (A.2), such that 
the expression in ( |5.1[ ) is fulfilled (approximately) at next-to leading order in perturbation theory. 
In particular, the next-to-leading corrections are reproduced within a relative accuracy of ~ 15% 
or better for all density- and velocity- correlations (and crossed components), independently of 
the linear input. Since the next-to-leading contribution is suppressed compared to the leading 
order, by a factor of a few in the relevant regime, this level of accuracy is also compatible with 
our estimated uncertainty in our truncated equation (of order ~ 10 %, as discussed in the text.) 
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